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Abstract —In a multi-user scenario where users belong to 
different operators, any interference mitigation method needs 
unavoidably some degree of cooperation among service providers. 
In this paper we propose a cooperation strategy based on the 
exchange of mutual Interference among operators, rather than 
of decoded data, to let every operator to recover an augmented 
degree of diversity either for channel estimation and multi-user 
detection. In xDSL scenario where multiple operators share the 
same cable binder the interference-cooperation (IC) approach 
outperforms data-exchange methods and preserves to certain 
degree the privacy of the users as signals can be tailored to 
prevent each operator to infer parameters (channel and data) of 
the users from the other operators. 

The IC method is based on Expectation Maximization esti¬ 
mation shaped to account for the degree of information that 
each operator can exchange with the others during the two 
steps of multi-user channel estimation and multi-user detection. 
Convergence of IC is guaranteed into few iterations and it does 
not depend on the structure of the interference. IC performance 
attains those of centralized receivers (l.e., one fusion-center that 
collects all the received signals from all the users/operators), with 
some loss when in heavily Interfered multi-user channel such 
as in twisted-pair communications allocated beyond 50-I00MHz 
spectrum. 

Index Terms —Digital subscriber line, multi-user receiver, EM 
iterative algorithm. Interference cooperation, distributed inter¬ 
ference cancellation, local loop unbundling, G. Fast. 

I. Introduction 

I NTERFERENCE mitigation is a largely investigated topic 
in communication systems for wired and wireless technol¬ 
ogy. The increasing demand of data-intensive services in fixed 
and mobile communication brings new challenges for inter¬ 
ference cancellation, and service providers (SPs) are looking 
for solutions that not only meet current, but mostly future 
requirements, still preserving the downward compatibility. 

Existing copper-wire infrastructure to provide net- 
bidirectional data rate of up to 200Mbps can use short 
loop lengths and high frequencies as in ITU-T G.993.2 
recommendation ID. Single cable binder can contain up to 
hundreds twisted pairs and these can be shared by multiple 
coexisting SPs. In multi-pair cable, crosstalk is the dominant 
impairment caused by the capacitive and inductive coupling 
among the twisted pairs. Techniques have been proposed 
to mitigate near-end crosstalk (NEXT) such as spectral 


shaping and frequency division duplexing El, 0. However, 
in digital subscriber line (DSL) the data rate is still limited 
by the far-end crosstalk (FEXT) that can be mitigated only 
by appropriate interference cancellation methods. More 
specifically, the performance is improved by using multi-user 
processing for FEXT cancellation, or vectoring 0. DSL 
Access Multiplexer (DSLAM) hosts signal processing units 
to mitigate upstream FEXT generated within a vectored 
group belonging to the same SR In common DSL scenario, 
non-vectored lines coexist along with the vectored groups in 
a same cable binder, and this prevents reliable FEXT control. 
Referring to multi-operator scenario in Figj^ vectoring 
technique is very efficient to cancel self-FEXT (i.e., FEXT 
of the same SP), but it has no control of FEXT due to the 
non-vectored lines and from other vectored group in the same 
cable binder (so called alien-FEXT) when twisted-pairs are 
shared among non-cooperating SPs. This problem is even 
more challenging in next generation G.Fast standard where 
the bandwidth is even larger than lOOMHz jSl over 50-200m 
cable length and alien-FEXT induced from just one temporary 
line-mismatch could loose all the benefits of vectoring at the 
expenses of energy 161. 

In this scenario, centralized vectoring where FEXT mitiga¬ 
tion is controlled by one single processing unit for multiple 
SPs would be beneficial but it is unfeasible due to the 
regulations of physical unbundling El. In addition, due to 
the privacy issue, SPs are not prone to exchange each other 
data to ease interference cancellation, even if in turn from 
this exchange there would be a benefit for all. Dynamic 
spectrum management methods could provide the same benefit 
as vectored transmission in the mix vectored and non-vectored 
(or equivalently non-cooperating SPs) cable environment 0. 
Noise decorrelation technique can handle alien-FEXT as ad¬ 
ditive structured noise, at the price of additional initialization 
(to estimate noise correlation) and adaptation (for time varying 
alien-FEXT scenario) algorithms El. If in vectored groups the 
non-vectored lines are unmanaged, the benefits of vectoring 
degrade rapidly Eol . The way to handle compatibility between 
vectored and non-vectored group lines is still an open issue 

El- 

Focus of this paper is the co-operation among SPs as 
envisioned in m to take the full advantage of vectoring 
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Figure 1. Interference Cooperation (IC) for multi-operator xDSL system. 


in multi-operator settings. Cooperation among multiple SPs 
(Figj^l being one alien-FEXT of the other is proposed here 
by avoiding the exchange of any sensitive data among SPs 
to guarantee that none of the other SPs can recover the 
information (either channel matrix and encoded data stream) of 
own users. The interference cooperation (IC) approach is based 
on the exchange of alien-FEXT only that is purposely removed 
as sketched in Eigj^ The iterative ICs that account either 
for multi-user channel estimation and multi-user detection 
are based on Expectation and Maximization (EM) estimation 
methods adapted to EEXT mitigation from cellular communi¬ 
cations miiniiii. In this paper, these iterative EM methods 
are reformulated by guaranteeing that every SP exchanges 
the interference with others after stripping the information 
depending on the own users to have an IC that preserves the 
sensible information. Given the linearity of the communication 
model, multiple-channel estimation and multi-user detection 
are handled using the same principle, with slightly differences 
due to the peculiarities of the two processing steps. Even if the 
signaling of the inter-operator interference seems to be a trivial 
extension of EM methods, the specihcities of interference 
cooperation need to address the overall multi-user channel 
estimation (comprehensive of alien-EEXT) and detection as 
a new problem. 

Benehts from cooperation in multiple receiver systems was 
established based on information theoretic capacity of multi¬ 
receiver cellular network 03 , ca, showing that the capacity 
loss due to inter-cell interference (i.e., alien-EEXT in wired 
systems) could be eliminated by cooperating base stations 
(BSs). In view to densify the BSs, multiple remote radio 
equipment are connected to a central BS through high speed 
link to get the benefit of central processing for up-link inter¬ 
cell interference mitigation ini. Interference alignment and 
cancellation has been proposed and validated in field trials 
ifTSl but still based on exchange of encoded data. Interference 
cancellation based on iterative subtraction of interfering sig¬ 
nals in cooperative BS clusters is in lfT9l . Il20l . Distributed BS 
co-operation scheme based on soft-combining for multi-user 
multi-cell system is in ED. In all these schemes, the degree of 
co-operation is controlled by limiting the backhaul capacity. 


Mixed soft and hard information exchange among multiple 
BSs for interference cancellation is in ll22]l . All these methods, 
in addition to E?l . E4l . ESll . Ebl . are excellent references 
to highlight the active research in the field of interference 
cancellation when multiple processing units need to cooperate 
by exchanging decoded (still sensitive) data to increase the 
spectrum efficiency in wireless and wired systems. However, 
this paper moves on a different conceptual setting where the 
exchanging among SPs is only the interference from received 
signals after iteratively stripping own sensitive data, or data 
mixed with an unknown mixing matrix acting as randomized 
multiplicative data perturbation ED- 

Security and privacy are primary concerns in communi¬ 
cation systems. Providing enough privacy should be also a 
challenge in any co-operative interference mitigation technique 
in multi-operator scenario due to the exchange of subscriber’s 
private information for the purpose of canceling the interfer¬ 
ence. Conventionally, security is provided at network layer for 
the well-known security threats e.g. eavesdropping, man-in- 
the-middle attack etc. ESll . E9]| . but the concept of security is 
changing in the advanced communication system design. Ex¬ 
amples of attacks on femtocells are in ESI. Recently, research 
community is giving more attention to physical layer security. 
Stochastic geometry approach for physical layer security in 
cellular system has been proposed in ED- This paper is 
motivated by the need to do all necessary to prevent the other 
SPs to decode the information on own users still having a 
beneht from the mutual SPs cooperation. 


A. Contribution 

Interference cooperation (IC) for multi-operator is the same 
for cellular and wired system, but in this paper we specialize 
the method for xDSL system, even if extension to multi-cell 
processing is conceptually straightforward with minor adapta¬ 
tions (not covered here). To better highlight the contribution, 
let us consider Eig[Dfor the simple case of two SPs with one 
user each labeled as 1 and 2 (channel estimation would be 
similar except for larger size systems, see Section III). Signal 
model is 

J/l = hiiXi+h2lX2+Wi=Sl^ + Wi 
2/2 = hi2Xi -I- h22X2 +W2 ^ g2^ + W2 


where SP 1 and 2 receive the interference from the other 
{h 2 iX 2 and hi 2 Xi), and wi,W 2 denote the additive white 
Gaussian noise (AWGN). In conventional data-exchange meth¬ 
ods each SP is aware (or estimated separately) of the channel 
vector gi = [hii,h 2 i]^ and g 2 = [hi 2 ,h 22 ]^, symbols xi 
and X 2 are locally decoded as xi and £2 (by SP 1 and 2, 
respectively) and iteratively exchanged to remove the crosstalk 
as yi — h 2 iX 2 to estimate xi by SPl, and similarly for SP2 
with j /2 - hi 2 Xi. 

In the proposed IC method, the overall system is modeled 


as: 



= hixi -I- 112 X 2 + w 


with hi = [/iii,/ii 2 ]^ and h 2 = [/i 2 i,h 22 ]^ the channels 
known to SPl and SP2, respectively. The interference is 
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exchanged at every iteration (i.e., based on the decoded xi, the 
residual j/i — huxi and hi 2 iiare sent 1—^2 that use by SP2 
for decoding £ 2 ', a similar reasoning holds true for signaling 
2 —> 1) to yield the 2 x 1 set of observations y — h 2 i 2 and 
y — hiii at SP 1 and 2 to decode ii and £2 using the locally 
known channels hi and h 2 . Differently from conventional 
iterative receivers where each exchange of data is to guarantee 
that the other processing units can reduce the interference, 
here each interference exchange favors the other SPs. Even 
if interference exchange seems less efficient than encoded 
data exchange, this paper proves that the iterative method 
solves implicitly the centralized problem of data detection 
(and similarly for channel estimation) by reorganizing the 
interference exchange to carry out Jacobi iterations distributed 
over SPs, that in turn converge to the centralized interference 
cancellation method within few iterations at price of signaling 
overhead. The convergence-rate of IC depends on the degree of 
interference that is faster when the channel matrix is diagonal 
dominant, but still within 3-5iterations. Even if the exchange 
among SPs is interference-based to preserve the privacy on 
own users, there could be a way for the other SPs to infer the 
information on the own users (e.g., SPl can extract h 2 i £2 after 
the stripping of SP2, but still mixed with unknown /121 and 
£ 2 ) by using blind separation methods Ea, 133]. However, 
blind-estimation methods need strict assumptions (e.g., known 
constellations and/or well algebraic structure of the channel) 
that make them unfeasible in practical xDSL systems (see 
Section III). 

Paper is organized as follows. Section II illustrates the 
system model for multi-operator scenarios with self-EEXT 
and alien-EEXT. Section III describes the co-operative channel 
estimation using iterative information shared among the SPs 
that implements the EM method distributed among different 
SPs with partial information on the interference. Section IV 
describes the cooperative multi-user detection (MUD) still us¬ 
ing iterative EM with MMSE estimator for multiple-operators. 
Numerical validations are in Section V, while Section VI draws 
some conclusions. 

II. System Model eor Multi-operator Receiver 

System model is for X-operators scenario with the channel 
matrices (or channels in short) for self- and alien-EEXT. 

A. Multi-operator System Model 

Multi-operator xDSL system for upstream is shown in Eig|^ 
where multiple service providers (SPs) sharing the same cable 
binder have N customer premises equipments (CPEs) labeled 
as CPEi i, CPE 2 ,i,..., CPEat i that are mutually interfering 
one another. To have an analytically tractable problem, here 
we assume that all the K SPs have the same number of 
CPEs mutually synchronized with the same upstream frame 
structure without any frequency drift one another so that 
any inter-carrier interference can be neglected, and we can 
employ a system model that accounts for EEXT on each carrier 
independently. 

Eig|^ shows the ensemble of all ITTV x 1 received signals 
y (f) at time t as given by 



Figure 2. System model for cooperative xDSL with K SPs and N CPEs 
each 


yi (t) 


r Hii 

H 21 

Hki ] 

r (t) 1 


r wi(t) ] 

y2 (*) 

= 
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+ 

W2 (t) 

yK (*) - 



f^2K 

Hkk . 
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- (t) . 


y(t) H x{t) w(t) 


where (i) = [xkp{t) ,Xk, 2 {l) ■ ■ ■ {t)]’' denotes the 

signals from the hnite alphabet A of M-QAM constellation 
that is transmitted from the N CPEs belonging to the fcth 
SP. The ensemble of NK x 1 transmitted signals related to 
all CPEs from all the K SPs are in x(t). Each N x N 
matrix Hii is the i —>■ j channel from N CPEs of the ith 
SP (CPEi i,... jCPE^r^i) towards the jth SP (/ — j) and it 
accounts for the channel including the self-EEXT (when i = j) 
and alien-EEXT (when i ^ j). H is the KN x KN channel 
matrix of the multi-user/multi-operator system. The AWGN 
Wfc(f) is uncorrelated among CPEs and SPs, with the same 
power: w (/) ^ CN{0,a'^T). Since the proposed IC method 
handles mutual interference as useful signals for the other SPs, 
for the scope of the paper the model ([T]) needs to be re-ordered 
by grouping the CPEs that belong to the same SP 

y (f) = HiXi (t) + H2X2 ...-y HkXk (/) + W (/) ( 2 ) 

where HI 2 , ■ ■ ■, ^ is the compound 

channel from the N CPEs belonging to the fcth SP 
(CPEi fc,..., CPEat fc) toward all the SPs as sketched in Pig|^ 
The arrangement © will be used by the fcth SP to estimate the 
KN X N channel in IC channel estimation (Section III), 
or in detection assuming is known (Section IV). Model for 
interference co-operation among the K SPs is shown in Pig|^ 
where all SPs share alien-EEXT interference information by 
using backhaul for interference exchange that let every SP 
forward to all the K — 1 SPs the received alien-interference 
after stripping the own data. Eor AT = 2 it reduces to a single¬ 
link between SPl and SP2. 

To separate channel estimation from decoding beyond the 
logical separation, the transmission is organized into frames 
with a frame structure that alternates a set of training samples 
with data to be decoded. The fcth SP assigns to its N CPEs 
T samples of the training sequences {x^ selected 

independently by the fcth SP. Training used by the fcth SP are 
not known to the other SPs to prevent estimation of alien- 
EEXT by other SPs. To ease the attention on IC method. 
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training can be considered as ideal (i.e., mutually orthogonal) 
even if some degradation is expected when this condition is not 
met. Needless to say that the channel estimation step includes 
the estimate of any frequency drift and their correction for the 
own CPEs, if necessary (not covered here). Channel estimation 
and MUD are discussed below as two distinct processing steps, 
still sharing several commonalities due to the linear model Q. 

B. Self and Alien channel models 

Channel model depends on the cable length and frequency 
range as increasing the frequency strengthen the inter-cable 
coupling, and this makes the FEXT comparable with losses of 
direct links (insertion loss). A simple model adopted here for 
numerical validations is based on the assumption that the direct 
links are normalized to have unit-amplitude, and the degree 
of diagonal dominance for self-FEXT channel Hu and alien- 
FEXT channel Hy (from {CPE^ ^ }^=i to Jth SP) is accounted 
in terms of scaling term a, the entries are 

with 2tt) (3) 

--CA/'(0,a^),Vp 7 ^ g (4) 

with the same coupling a for both self- and alien-FEXT. 

Off-diagonal terms are as small as a < —IQdB up to 
10-30MHZ bandwidth 01 Ql, these values are compliant 
with channel measurements of 24 pairs (0.5mm/ea) over 
200m twisted pair cable up to 80-100MHz bandwidth as for 
next generation DSL ll35]l . The FEXT coupling a increases 
vs frequency and for long-cables, it is not so unusual to 
have a FEXT that is comparable to direct link when above 
50-lOOMHz as for G.fast settings O) adopted in numerical 
validations (Section V). 

III. Cooperative Channel Estimation 
The training sequences are composed of a set of T complex¬ 
valued samples (T > N) transmitted by all the KN CPEs at 
the same time. The N x T training sequence matrix for the 
N CPEs of fcth SP {CPE^^fclfci is 

Xfe = [xfc(l),Xfe(2),.-.,Xfe(r)], (5) 

over-bar denotes the frame of training sequences. These train¬ 
ing are pseudo-random sequences assigned by SP to every 
CPE and never exchanged among the SPs so the A:th SP 
knows Xfe, but not the others. To simplify, it is assumed that 
the training sequences of one SP are mutually uncorrelated 

with training sequences of other K — 1 SPs and have ideal 
_ _ ^ 

auto-correlation so that Xi;Xj, = I and X^Xj = 0 for 
all i k. Notice that for random selection of the training 
sequence as in practice happens, these ideal properties hold 
only approximately at price of a negligible loss of perfor¬ 
mance (Section V). The discrete signal received at fcth SP 
Yfe = [yfc(l),yfc( 2 ),...,yfc(T)] is 

Yfc = HfcfeXfc + HmfeXm -b Wfe (6) 

m^k 

here separated into own data (Hfe/jX^) and alien-FEXT 
( HmfcXm). Below the conventional Data Cooperation 

m^k 

(DC) is discussed first as reference scenario, and then the 
method based on Interference Cooperation (IC). 


A. Data Cooperation (DC) 

In DC the initial signaling setup lets all SPs be aware of 
the training sequence of all the other SPs. In this case the 
channel estimation can be carried out by each SP to attain 
the centralized approach where virtually one SP collects all 
the received signals. More specifically, when the /cth SP is 


aware of all training X = 
becomes: 

Yfe = G 


•- X - X 


feX 


2 5 • 


■Wu 


—T 

,x^ 


n T 


the model 


(7) 


where = [Hifc, • • • , Hkk, ■ ■ ■, collects the channels 

from all CPE^ ^ toward fcth SP. Maximum likelihood estima¬ 
tion (MLE) of this compound channel is 

Gfc = YfeX"(XX^)-i (8) 


According to the properties of training sequence, the covari¬ 
ance of the channel estimate G^ is 

2 

cov(Gfc) > '^I, (9) 

that is the Cramer Rao bound of the channel estimate. These 
estimates are independent and error scales with the training 
sequence length T and AWGN power cr^. Once again, it is 
crucial to remark that DC mimics the centralized approach 
so that fcth SP estimates self-FEXT (Hfcfc) and all the alien- 
FEXT (Hmfc for Vto 7 ^ fc) channels that are locally used by 
fcth SP in MUD. If channel is slowly varying as in wired 
system, this method needs a minimal inter-SP signaling at 
setup to have a consensus among the usage of training and 
the channel estimations can be carried out independently by 
each SP according to ([^. 


B. Interference Cooperation (IC) 

IC method is iterative and implies not to exchange the 
training among the SPs (as considered a sensitive information) 
but rather the training-induced FEXT. The model for channel 
estimation by the fcth SP as the SP of interest is by collecting 
the T training samples into Y = [y(l), y(2),..., y(T)] 
according to ( 0 : 

Y = HfcXfe + ^ H„X^ + W, (10) 

m^k 

where contribution by the different CPEs are grouped ac¬ 
cording to their SPs. Since each SP is only aware of its 
own training sequence, after the fcth SP receives the alien- 
FEXT signaling from all the other K — I SPs in the form of 
Ym — X HimX/ (say from the mth SP withm 7 b fc), it stacks 

l^k 

all the received alien-FEXT to locally reproduce the following 
model complementary to ( [T 0 | ) 

Y - ^ H^X™ = HfeXfc + W, (11) 

m^k 

after the signals from all interfering SPs are stripped out from 
each SP and forwarded to the fcth one. The linear model •HD 
is solved iteratively in the IC method as detailed below. 

Let be the initial estimate at fcth SP by assuming that 
alien-FEXT is an augmented AWGN, the fcth SP exchanges 
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Yfe — with the other K — 1 SPs in trade of — 

HmmXm one by one from all the others SPs (here m ^ k 
denotes any arbitrary SP different from the fcth one). At the 
nth iteration of the iterative IC method, the fcth SP estimates 
re-encoded as H^^X^ and forwarded to 
the mth SP in form of alien-FEXT, and similarly each of the 
other K — 1 SPs re-encodes H^^X^ for the beneht of /cth 
SP The interference is thus exchanged among all the K SPs 
one by one as N x T matrices, kth SP uses the received re¬ 
encoded interference for alien-FEXT cancellation and iterative 
channel estimation based on the linear model GD adapted 
for the iterations. It must be noticed that at every iteration 
the estimate by the fcth SP is the channel KN x N matrix 
Hfc = [Hji, H^ 2 ) • • ■) ’^Ik] ^ stacks not only self-FEXT 
but also alien-FEXT channels from own CPEs towards 
the other SPs for m ^ k). At every iteration, the cost 

function is 


H 


("+!) 


= arg min 

Hr, 


ry{n) iqr -y 


( 12 ) 


where = Y— ^ Hm^X,„ is the effective variable that 

ra^k 

collects all interference Hm^X^ exchanged by the K—\ SPs 
(m ^ k) locally used for rehnement of H/j. The estimation of 
Hfc becomes 

h(.+i) , z<">xf (XiXf )-■ (13) 


Once again, SPs cooperation rehnes self-FEXT and alien- 
FEXT by constraining the cooperation among SPs to never 
exchange the training X^ and X^ among each other, at 
most mixed by the (unknown) channel responses. The iterative 
channel estimation algorithm is summarized in Algorithm-1 
for the kth SP. 


Algorithm 1 Upstream Multi-operator Interference Coopera- 
tion Channel Estimation_ 

. Initialize =YfeXf(XfeXfri 

• Receive alien-FEXT from all K — 1 SPs 

(m 7 ^ k). 

• Interference exchange (m <->■ k): evaluate = (Y^ — 
H^LX™)xf (XfeXf and forward X^ to mth 
SP in exchange of H^^X^. 

• for 72 = 1 : Niteration 

- Estimate 

H<r‘> = IY. -h2x„ I xf (x,xf )-■ 

y m^k j 

- Receive alien-FEXT Y^ — — 

^ HpmXp from all AT — 1 SPs (m ^ k). 

pz^m,k 

- Interference exchange (m o k): evaluate = 

(Y„-H^+1)X„- ^ H(’^Xp)xf(XfcXf)-i 

p^m,k 

and forward to mth SP in exchange of 


1) Convergence of IC: Iterative IC method for channel 
estimation is equivalent to iteratively solve the linear system of 
Muations for the (centralized) model Y = HX -t- W from the 
(hjl. In other word, given the centralized model Y = HX-pW, 
the training can be partitioned into self and alien-FEXT 
training (X = Xg-I-Xa) so that the IC can be rewritten as Y = 
H("+i)Xs -PH^^^Xq -PW, this is the basis to prove that the 
iterations in IC is equivalent to Jacobi iteration 13^ for channel 
estimation such that the estimate is carried out independently 
by each SP without exchanging X^ but rather in form of 
mixed values H^"^Xa. Some technicalities are necessary for 
the equivalence. Let M = blockdiag [Xi, X 2 ,..., Xk] be the 
block diagonal of the training sequences known by each SP, 

—T 

where Xk = Ik ® Ia^ here reordered just to comply 

with matrix algebra computations, the block off-diagonal term 
N = M — (\k ® [<^ 1 ) <^21 ■ • • j ^k\) denotes the matrix with 
training for alien-FEXT, the IC method (Algorithm 1) reduces 
to the set of Jacobi iterations 

M • vec = Nmec -P (l^r O uec(Y)) (14) 

distributed over the SPs as each SP iteration rehnes (locally) a 
portion of the channel estimate vec where the estimates 

by each SP are arranged into a vector by vec{.) operator. Proof 
of convergence of Jacobi iterations is in Appendix-A, but since 
Jacobi iteration converge to = YX^(XX^)“^ for n 

large enough (in practice, after 3-5 iterations, see Section V) 
the IC converge to the centralized MLE of the channel matrix 
without any exchange of data/training except as mixed values. 
It is crucial to remark again that in IC method the channels 
Hfc = , H ^2 j • ■ ■) ^ that cause alien-FEXT towards 

the other remaining AT — 1 SPs is the one that is updated within 
the kth SP. At convergence, this information is resident in fcth 
SP in form of channel estimates, and this enables the IC Multi¬ 
user Detection (IC-MUD) as detailed below. 

IV. Cooperative Multi-user Detection 

Data detection in multi-operator environment is based on 
channel state information available at each SP and the degree 
of cooperation among SPs. To ease the analysis still consistent 
with the problem at hand, the channels are assumed as random 
and mutually orthogonal i.e. E = 0,Vi f j. The 

signals received at fcth SP for data detection follows from Q 
(time index t is omitted to simplify the notation): 

Yk — ^^kk^k H" ^ ^ '^mk^m ^k (15) 

m^k 

where HfefcXfe, Hmk^m and accounts self-FEXT, alien- 
FEXT and AWGN. The iterative MUD for Data Cooperation 
(DC) and Interference Cooperation (IC) are based on the EM 
method El, El, ia, ED, and these iterative methods are 
compared to the centralized MUD (i.e., all the SPs are decoded 
jointly) used as a reference scenario. 

Once defined the linear model, MUD can be based either on 
zero-forcing (ZF) or minimum mean square (MMSE) criteria 
depending on the diagonal dominance of the channel model 
El . Herein all MUD methods are based on the matrix decision- 
feedback equalizer (DEE) from QR decomposition of the 
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corresponding channel (or equivalently, BLAST system) either 
using ZF or MMSE criteria as widely investigated and adopted 
in vectoring (see e.g., |4^ and ID). To simplify the notation 
of MUD, the metrics are referred to ZF-MUD and matrix DFE 
are left indicated as well established. 


A. Centralized Multi-user Detection 

In centralized MUD all SPs forward the received signals 
to a central processing unit that is aware of all channels H 
(possibly estimated) and jointly decodes all the data streams 
according to the system model ([T]). The ZE-MUD is 

X = H-iy, (16) 

and the variance of the noise at decision variable is lower 
bounded by <7^/(1 + {NK — l)a^), it decreases with the 
number of cooperating SPs {K) and CPEs {N) due to the aug¬ 
mentation of signal flows in vectoring. For larger bandwidth 
(say >50MHz) the FEXT-coupling increases and the MMSE- 
MUD with matrix-DEE takes more efficiently into account the 
higher FEXT in vectoring. Signaling constraints for data-rate 
and privacy issues make the full cooperation with centralized 
MUD unfeasible but still it is an excellent reference bound in 
performance analysis (Section V). 


B. Data Cooperation (DC) MUD 

In DC iterative technique, it is assumed that the channel 
Gl = [Hife,...,Hfcfc,...,Hiffc] is known to fcth SP, or 
estimated as a separated step discussed in Section III-A. At 
every iteration (say nth), every SP forwards its decoded data 
Xm^ to all the other K — 1 SPs that in turn cancel locally 
the alien-EEXT using all the received data from the K — 1 
SPs after channel reshaping this step employs the 

following model for every iteration 

Yk- = H/cfcXfe + Wfc ( 17 ) 

m^k 


The ZE-MUD simplifies into a set of soft-detection for each 
decision variable as ini: 


(n) 

H =9a 



(18) 


where pA[a;] = Ex^a{x\x -f z} is the symbol-based MMSE 
estimator at the decision variable and is the correspond¬ 
ing (iteration varying) variance of the noise at the decision 
variable. The soft-decision pa[.] is modulation-dependent as 
described in Appendix C, and it reduces to tanh[.] for BPSK 
and QPSK modulations Cl. Smoothness depends on 1 /ct^ 
and it can be “hardened” (by reducing ct^) at last iterations, if 
necessary, or by updating its value as in The same soft- 
decisions are employed within the loop of the matrix-DFE to 
avoid error propagation within the same SP. 

At 1st iteration, data from other SPs are not available at any 
SP, therefore from the received signal the noise at decision 
variable is assumed Gaussian and alien-FEXT augments the 
AWGN: ^ + w^) ^ CN (O, afl) with a? = 

m^k 


((iT — 1) Na'^ -f CT^) /(1-|-(A^—l)a^). Once the decoded data 
is shared among the SPs, the variance at decision variable 
progressively reduces down to the (almost) complete alien- 
EEXT cancellation cr^ = <7^/(1 -I- (A^ — 1 )q;^), for n large 
enough. However, it must be noticed that the noise variance 
in DC is always higher than the centralized approach as in DC 
only the N self-FEXT links are used for data estimation rather 
than in centralized vectoring approach where all KN links 
contributes in decoding. Compared to centralized vectoring, 
the degradation of noise power at decision variable for DC is 
at least 

N{K-l)a^ 


DC 


loss 


= 1 + 


1 + {N-I)a'^ 


(19) 


that for a very large number of CPEs it depends on K as 

DC/oss — K. 


C. Interference Cooperation (IC) MUD 

For IC method, the reference system is (|^ and the ensemble 
of channels HI 2 , • ■ •, is known to fcth 

SP as being estimated by the IC channel estimation (Section 
III-B). At each iteration, all the K SPs estimates the sym¬ 
bols x^”\ X 2 ”\ ..., x^^ and forward each-others the received 
signals after stripping the own signals and leaving the alien- 
FEXT. Eocusing to the /cth SP, it stacks all the measurements 
received by all the K-1 SP into the KN x 1 vector 

= (20) 

m^k 

that is the results of stripping Xm^ with the local KN x N 
mixing channelHm known only to the mth SP (m k) as 
detailed in Section III. The model at the fcth SP reduces to 

= HfcXfe-f w, (21) 


that is made redundant hy KN lines as the result of the 
augmentation from the alien-EEXT exchange. 

Signaling exchange and MUD processing for © is detailed 
in Algorithm 2 for ZE criteria (MMSE differs only for the 
metric), and it is illustrated by focusing to a pair of SPs, say 
TOth and kth. The mth SP estimate Xm^ is mixed with (locally 
available) channel to have that is forwarded to 

the kth SP, and symmetrically from fcth SP so that the inter¬ 
ference is exchanged one-by-one among all the K SPs. Based 
on these exchanges, fcth SP mitigates the alien-EEXT from the 
locally received signal y^ as y^ — ^ Similarly, 


m^k 

the TOth SP estimates the alien-EEXT contribution from fcth 

using own CPEs (HmmXm^) and 


SP as y„ 


- X) Hi^xJ 

l^k 

alien-EEXT (Hi^X;"^ for I k,m) received from the other 
SPs. The estimated alien-FEXT is again exchanged among all 
the K — 1 SPs and the fcth SP to locally build the new KN xl 
vector (201 for MUD x^”'*’^^. 

After mutual alien-EEXT signaling, the estimation of 


(n+l) 


from the linear model (21 1 is carried as part of the 


EM iterations na, El. Namely, the cost function for ZF- 
MUD (similarly for MMSE-MUD) is given by 


X 


(n+l) 

k 


= arg min 

Xfc 





2 


( 22 ) 
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and it is be solved by matrix-DFE from QR decomposition of 
the compound KN x N mixing channelHfc (see Algorithm 
2 for details). Due to the iterative exchange of alien-FEXT 
in interference cancellation, it is convenient to adopt the soft- 
estimates to avoid error propagation and have a gain in SNR 
of approx. 5dB (Section V) in line with ref. ll24l . The MUD 
of the symbols at (n + l)th iteration becomes; 


(n+l) 


= 5A 




(23) 


where tr^ is the variance of the noise at decision variable 
for IC method and pa[-] is symbol-based MMSE estimator 
based on the alphabet A (Appendix C). Variance cr^ reduces 
vs iterations and it can be calculated while evaluating the new 
soft-decisions El. 

At the setup, there is no information available on alien- 
FEXT and thus the estimate (23 i is replaced by the estimate 
based on yfc that accounts for alien-FEXT as an augmented 
AWGN. After the 1st iteration, IC and DC differ as received 
signals are stripped of the own data to form the exchange 
signal Yk — as detailed above (or in Algorithm 


m^k 

2). After alien-FEXT exchange, the noise power at decision 
variable is bounded by cr ^/{1 + {NK—I)a'^) as for centralized 
approach. 


Algorithm 2 Upstream Multi-operator Cooperative ZF-MUD 

• QR decomposition of Hkk (iteration 0): QR = and 

Yk = Q^Yk- 

• Initialize 



[yfc]i 


N 




j=i+l 



, for i = A^, A^ - 1,..., 1 

• QR decomposition of (iteration n > 1);QR = Hfc, 

• for 77. = 1 . Niteration 

- Receive from allAT — 1 SPs (to ^ k). 

- Interference exchange (k o to): evaluate Yk — 

^ and forward to TOth SP in exchange of 

l^m 

Ym- H/mx[”^. 

l^k 

- Estimate 


(n+l) 


5A 




N 

7 ^ . ■ 

- -11 ' 
]=i+l 


(n-l-l) 


J 3 


, for z[,"^ = 

- Update based on 


Fig[^ illustrates the (n + l)th iteration of IC MUD Al¬ 
gorithm 2 for K — 2, there are two signaling phases with 
an overall signaling complexity Os (2A^), and computational 
complexity Oc (4A^^) at each SP. Generalization to arbitrary 
K SPs makes the overall complexity Os (2 {K — 1) N) and 



Figure 3. Signaling and Computational Complexity at each SP for K = 2. 


Oc [2KN‘^y Compared to DC MUD with signaling complex¬ 
ity Os {{K — 1) N), the signaling of IC MUD is twice with 
remarkable benefits (Section V). Even if the data-rate is quite 
intensive, the SPs are hosted in the same cabinet on at short 
distance, and there are several parallelism that can be exploited 
in this setting for any practical implementations M- 

1) Convergence of ZF-MUD IC: Similarly to channel esti¬ 

mation (as IC MUD rely on the same linear model Q used in 
multi-user channel estimation), the updating ( |2^ is conceptu¬ 
ally equivalent to the partition of the channel model into self 
(HQ and alien-FEXT (HJ so that the iterations reflect this 
partitioning eased by the IC signaling, and the local availability 
of channel matrices Hfc: y = -f -f w. 

The proof of convergence toward the centralized estimate 
follows the same steps that highlights the equivalence between 
iterations •HD and the Jacobi method to solve the linear 
system y = Hx-fw |[36l. By reshaping the linear problem 
into block diagonal matrix M = blockdiag [Hi, H 2 ,..., Hr] 
that accounts for self-FEXT (say H^), and block off-diagonal 
matrix N = M — {Ir ® [Hi, H 2 ,..., H^-]) for alien-FEXT 
(say Ha), each IC MUD is equivalent to the Jacobi iteration: 

Mx(”+i) = Nx(") -f {\r ® y) (24) 

carried out separately by each SP as MUD step (j2^, except for 
soft-detector pA [•] that is unavoidably related to the alphabet A 
and used to prevent error propagation. Once again, the proof 
of convergence of Jacobi iterations (Appendix B) guarantees 
that IC converges to the same centralized vectoring solution 
for ZF metric. It can be shown that proof for MMSE-MUD 
IC follows the same steps. 

2) Privacy: The exchange between any two SP is mixed 
by alien-FEXT channel that is never know according to the 
IC method for channel estimation (Section III). The fcth SP 
could recover the mixed signal HmfcXm^ but this is without 
the knowledge of the mixing alien-FEXT channel H^fc that 
in IC method is not known to the fcth SP, but only to the 
TOth SP, and this randomizes Xm to preserve its the privacy 
Ea. Even if the exchange among SPs is interference-based, 
in principle there could be a way for the other SPs to infer 
the information on the own users. For the example at hand, 
the estimation of the alien information Xm could reduce to a 
multi-user blind separation from H^fcXm^ with the unknown 
mixing matrix Timk- However, even if the modulating terms 
are non-Gaussian as a necessary condition for blind-estimation 
methods, these are unknown in term of constellation, power 
(or scaling factor) and phase-stationarity (e.g., each user can 
have an arbitrarily rotated constellation specifically employed 
to prevent the decoding by alien SPs, but still being determisti- 
cally known for the own SP) to prevent a reliable decoding of 
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N X N MIMO mixing that, for the alien-FEXT channel Hmk, 
it is likely below the decoding capability. In this sense the 
privacy can be considered pragmatically preserved to enable 
the adoption of IC methods from commercial operators and 
enable local loop unbundling. 


V. Numerical results 

Numerical simulations validate the asymptotic performance 
of the IC method and evaluate the effectiveness of the coop¬ 
erative method based on the exchange of interference only. 
We consider a scenario with K = 2 and 3 operators with 
= 10 CPEs each as being a reference for G.fast setting. 
Performances of the IC algorithms are evaluated either for 
statistical model ifJTl by varying the FEXT coupling a, and 
by considering the cable-models of 100m length ll35l up to 
the frequency of 200MHz to simulate the behavior in G.fast 
settings. In any case, the direct link is independent over the 
lines (all lines with the same length) so that HkklPjP] — 
exp(j0p) with dp ^ U{0,2tt) statistically independent over 
lines. The transmitted symbols from every CPE belong to 
QAM constellation with a size that depends on the degree 
of G.fast specifications. Power spectral density of the signal at 
CPE is -76dBm/Hz and the noise is -140dBm/Hz as customary 
irSl . Il38]l . 091 . signal to noise ratio (SNR) is always referred 
at the receiver unless defined at decision variable. 


A. Iterative Channel Estimation 


Training are generated from a random set of QPSK symbols 
known to each SP as this approximates the ideal training with 
XfeXf ~ I and XfcXf ~ 0 for all i ^ k. PigQ shows 
MSE vs iterations of self (H^Q and alien-EEXT (H^) channel 
estimates for K = 2 and 3 operators scenario using IC, and 
the corresponding Cramer Rao Bound (CRB) (j^. The perfor¬ 
mance has been evaluated for SNR= G {5dB,15dB} 
and EEXT coupling a G {0.5,1.2} to validate the con¬ 
vergence of the proposed algorithms in these severe EEXT- 
conditions. The MSE is normalized by the power of channel 


for self-EEXT E 


IH,, 


> E 


IIH, 


and alien-EEXT E 


|H, 


since 


the normalized MSE for Hi 


lower than MSE of H,i. Centralized method attains the CRB 


MSEcrb = trace 


•-IJ 

a2(Xx")-i 


and here IC attains the 


CRB in 3-5 iterations for K = 2. When K > 2 a larger 
number of iterations is necessary as the iterations are implicitly 
pairing the alien-FEXT to each users during the alien-EEXT 
exchange steps. If EEXT-coupling a is smaller, the number 
of iterations for IC convergence decreases as expected for 
smaller interference settings. In G.fast scenario the IC method 
converges (with in H-2-3dB of excess MSE compared to CRB) 
in 2 -t 3 iterations up to a frequency of 100-120MHz, even if a 
practical convergence (not necessarily to the CRB) within 1-2 
iterations. The cost of inter-SP signaling in this case is limited 
to the set-up phase and it needs not further optimizations to 
comply with inter-SP data-rate. 



a. [dB] 


Eigure 5. SNR at decision variable of MMSE-MUD methods vs a for 
K = 2 (empty marker) and K = 3 (filled marker) with SNR = 40dB and 
6 iterations. 


B. Iterative Detection 

Performance for multi-user/multi-operator detection are 
evaluated in term of SNR at decision variable (SNRu) for 
different TEXT coupling and modulations. MMSE criteria 
for MUD is employed to cope with large EEXT coupling. 
Eig|^ shows the SNRd for varying a and SPs AT = 2,3 
for SNR = 1/cr^ = 15dB scenario after 6-iterations for 
DC and IC methods with QPSK modulation. Performance 
of IC MMSE-MUD is compared with DC MMSE-MUD and 
MMSE-MUD with local vectoring (without alien-FEXT com¬ 
pensation) thus showing that the IC MMSE-MUD outperforms 
all the methods. Degradation with respect to centralized MUD 
is mainly due to error propagation for large EEXT coupling, 
soft-decoding in IC MUD prevents the propagation of errors 
and guarantees SNRd > lOdB even for a — OdB and 
K = 2,3. Any inter-SP cooperation becomes mandatory for 
large degree of coupling, say a > —20dB, with IC MUD 
uniformly when the error propagation dominates. Since in 
G.fast the crosstalk can be considered as a > —20dB, the 
cooperation among SP should be considered as mandatory to 
make an efficient usage in any condition of coexistence of 
users within the same cable bundling. 

Convergence analysis of IC MMSE-MUD for QPSK with 
soft (empty-dots) and hard (filled dots) decisions DEE is in 
Eig|^ showing the SNRd vs iterations for K = 2,3 and 
varying a. Once again, this numerical results validates that 
the iterative IC algorithm converges into few iterations (say 
3-5 iterations) for most practical purposes, with a remarkable 
benefit when using soft-decisions ( |2^ . Convergence is faster 
for smaller a and slower for large number of SPs (here K=3) 
as iterations implicitly ease the assignment to each SP the 
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Figure 4. MSE vs Iterations(n) of IC channel estimation in multi-operator xDSL for SNR = [5,15]dB,a = [0.5,1.2], K = 2 (left) and K = 3 (right). 



123456789 10 

Iteration 


Figure 6. SNR at decision variable vs iterations of IC MMSE-MUD for 
K = 2 (empty marker) and K = 3 (filled marker), DEE with soft (dashed 
line) and hard (solid line) decisions, and varying a =(-10(circle),0(square)) 
dB. {N = 10, SNR = IbdB) 


alien-FEXT corresponding to the own users to be exploited 
efficiently as “useful signal”. 

Average performance from all lines vs frequency for G.fast 
channel model Ea, 123 (cable length is 100m for all the 
CPEs, N=10 CPEs per SP and K=2 SPs) is in Pig|^ channel 
model is sketched on top-corner for insertion-loss (IE) and 
EEXT, and the corresponding a. The thresholds on SNRo 
to attain the symbol error rate of 10“^ are superimposed for 
different QAM constellations, from BPSK to 4096-QAM. The 
MMSE-MUD for centralized system is the upper-bound with 


some dispersion of values (shaded area) but still above the 
BPSK threshold. The lack of any alien-EEXT mitigation as 
for no-cooperation scenario offers the lower bound. The DC 
MMSE-MUD method shows comparable performance with 
centralized MUD up to approx. 120MHz (or equivalently from 
the channel parameters a < —IQdB) in accordance with 
threshold effect in Eig|^ IC MMSE-MUD (4 iterations) with 
soft-decoding outperforms the DC MUD and the same IC 
MUD with hard decoding for QPSK and 16-QAM (other 
modulations are consistent with these results), thus confirming 
all the results in Eig|^ In addition, the IC method guarantees 
the transport over all the G.fast bandwidth up to 212MHz. 

The SNR lets to compute the average throughput (Table 1) 
according to the “gap-formula” for bit-loading usually adopted 
in DSL lim, ||39l for SNRd at ^th tone 


= log2 1 + 


SNRdi 


where the gap U=10.8dB=6dB(SNR margin)H-9.8dB(SNR 
gap)-5dB(coding gain) according to the ITU specifications for 
symbol error probability 10“^, maximum loading of 12bits 
and framing overhead is 12% ll38]l . Over the bandwidth 2- 
106MHz DC and IC attain the same throughput of centralized 
MUD, but when extending the bandwidth up to 212MHz 
the loss of IC MUD is negligible compared to centralized 
method, and the DC looses approx. lOOMHz in throughput. 
The same Table shows that is it remarkably better to partition 
the overall bandwidth among the 2 SPs (e.g., alternate usage 
of the tones per SP) rather than allocate all the CPEs over the 
same bandwidth without any degree of cooperation, in this 
case the throughput is half the value for centralized (or IC) 
MUD. 


VI. Conclusions 

Conventional interference mitigation techniques are based 
on centralized iterative interference cancellation that exchange 
data without any limitation on ownership. In this paper, we 
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Table I 

Average throughput of different cooperation schemes for G. Fast settings (N=10, K=2, cable: 100m) 


Cooperation Scheme 

2-106 MHz 

2-212 MHz 

Centralized 

718 Mbps 

884 Mbps 

Equally shared bandwidth 

357 Mbps 

425 Mbps 

Interference Cooperation (IC) 

714 Mbps 

854 Mbps 

Data Cooperation (DC) 

711 Mbps 

756 Mbps 

No Cooperation (NC) 

275 Mbps 

275 Mbps 



Frequency, MHz 

Figure 7. SNR at decision variable vs frequency in G.fast settings (cable 
length is 100m for all the CPFs, N=10 CPFs per SP and K=2 SPs) for MMSF- 
MUD: IC method with soft (SD) and hard (FID) decoding, centralized MUD, 
DC method. Channel losses are in top-right comer. 


propose multi-user receivers employing channel estimation 
and detection in the multi-operator xDSL unbundling scenario 
where cooperation among different operators is obtained by 
exchanging the alien-interference considered as useful signal 
by the other operators. The IC method is iterative as interfer¬ 
ence exchange refines alien-interference and own data, and it 
has the merit that this cooperation among different service 
providers never exchange the decoded symbols considered 
as sensitive-information. Since each operator receives the 
interference from the others in exchange of their own, it can be 
considered an interference-based cooperation that can exploit 
all the redundancy from all the data-paths with respect to 
the other operators and thus the performance attains the same 
performance as a centralized vectoring method. 

The convergence of IC is impaired by error propagation and 
soft-decision multi-user detection guarantees the convergence 
in few iterations (type.3-5 iterations) with a threshold in loss 
when overall interference becomes comparable with signal 
of interest. Even if data-rate exchange among operators is 
double the interference mitigation method by data exchange, 
the parallelism over frequency can be easily exploited in inter¬ 
operator communication using specific architectures ll43l . 
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Appendix A 

IC Channel Estimation Convergence 

The exchanging of the interference reduces the channel 
estimation for K SPs to the linear system (except AWGN that 


is irrelevant for MLE method) Y = HX, after vectorization 
this is rewritten as linear system 


b = Ah 

to be solved with respect to h = vec (H) G C 


Terms in (25 1 are 


b = lif 0 vec (Y) 
as Xk = 


C 


(25) 

K^NT-kI 


where training 

Ik ® X^, C) 


A =\k ® [Xl^ X 2 

for the ^th SP are reordered 
In G n ^ x]^ of size K x 1. 

Block-matrix A of training sequences can be decomposed 
based on the knowledge of each SP as block-diagonal M = 
blockdiag [Xi, X 2 , ■ ■., Xk] (for each SP) and block off- 
diagonal matrix N = — (A — M) (signaled by the other SPs) 
so that the updates at every SP is equivalent to the Jacobi 
iterations: 

= Nh(”) -f b (26) 


to be applied in least-square sense as A is not a square 
matrix: h 

h("+i) = (M^M)”^M^Nh(”) -f b. (27) 

Convergence to = (A^A) ^ A^b regardless of 

the initialization depends on the spectral radius of 

(M^M) Since M is diagonal dominant, the 

spectral radius p((M^M) ^M^N) < 1, and then the Jacobi 
iterations always converges for any starting vector Proof 
is given below. 

Let is the error at nth iteration. As 

= Nh^"^ + b, the iterations are equivalent to the 
update and convergence is 

guarantee if the eigenvalues of (M-^M) are strictly 

smaller than 1. Rewriting the updating as 

(M^M)”^M"N = I-(M^M)”^M^A (28) 

where M^M = blockdiag [X(^Xi, X^X 2 ,..., X^Xk] , 



r xi^xi 

Xl^X2 ■ 

■ XI^Xk 1 

11 

ti: 

.. 

.. 

.. 


1 

X^X2 ■ 

1_ 


the update (M^M) ^M^A becomes: 

(M^M)”^M^A = 

I (xf X2 ■ {xi‘Xi)-'^xHxji 

{^XH X2)~'-xH X, I ... [xHx^Y'^xHXk 

- {xSxK)~^xgXi (xgxY~^xgx2 ■ I 

_PO' 

Since the training sequences are orthogonal (i.e., X^X^ = 
I and XfcXj = 0 for all i Y k), or at least uncorrelated, the 
entries of the off-diagonal blocks (^X^Xk) X^Xm are very 
small (ideally null). Hence the matrix is diagonal dominant and 
spectral radius is p((M^M) M^N) < 1, the convergence 
is guaranteed. 
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Appendix B 

IC Data Detection Convergence Conditions 

The exchanging of the interference reduces the detection for 
K SPs to the linear system y = Hx, after vectorization it is 


b = Ax 


(31) 


Terms in (31 1 are 


H e and 


to be solved with respect to x € 
b = 1^ 0 y e A = Ij 

Ik = [1,1, ■ • •, 1]^ of size A x 1. Block-matrix A of channel 
can be decomposed based on the knowledge of each SP as 
block-diagonal M = blockdiag [Hi, H 2 ,..., H^] and block 
off-diagonal matrix N = M — A so that the updates at every 
SP is equivalent to the Jacobi iterations: 

Mx(”+i) = Nx(") -f b (32) 


Similarly to Appendix-A, the iterative method is 

x("+i) = (M^M)”^M^Nx(") -f (M^M)”^M^b (33) 

System convergence to x*^°°) — (A^A) ^A^b that for the 
structure of the problem coincides with ZF MUD (18), and 
depends upon the spectral radius of (M^M) ^M^N. The 
structure of the matrix is: 


(M^M) = I - (M^M) ^M^A (34) 

where M^M = 6/ocfcdwp [Hf Hi, H^Ha,..., Hf 
and 



r HfHi 

HfH2 

HfH^ 


II 

H2"Hi 

H|^H2 

HI^Hk 

(35) 


HfHi 

HfH2 

Hf Hk . 


thus (M^M)' 

^M^A becomes: 




(m"m) ^M^A = 


I (HfHl) ^HfH2 ■■■ (HfHi) 

I ... 

( 36 ) 

Since the off-diagonal terms are random (3-5) and 

HfH™ = ~ 0 for large N and K, 

1=1 

(HfHfc)”^HfH^ is close to 0. Hence (M^M)”^M^N 
is diagonal dominant with p((M^M) ^M^N) < 1 and 
convergence is guaranteed. 


Appendix C 
MMSE Estimator gA[-] 

Soft-detector pa[-] plays a role to avoid error propagation 
and it depends on the transmitted constellation A. During 
iterations, each symbol for each user at decision variable can 
be modeled y = x + z: the sum of a complex valued symbol 
X G A and a Gaussian noise z ~ CA(0, cr^) with power 
that depends on the iteration. The conditional expectation 


g[y] — E{x\y} depends on the probability density function 
(pdf) py \y\ of y, and in turn on pdf of x and 2 ; as being both 
random variable statistically independent 

5A = v + [ 2 ^]} ■ 

Eor separable rectangular constellation (e.g., M-QAM: A = 
Afi X A/ with Afi = A/ = {zbl,ib3,..,ih(2'\/ M — 1)}) the 
soft-detector is separable pa[2/] = ((’[ur] + j’Pbji] onto real 
and imaginary component, and it resembles a soft-detector 
for multilevel constellations and it becomes hard-detector for 
-G 0 . 

To simplify, let x G {ai, 0 : 2 , cim} where M sym¬ 
bols of the alphabet are equally likely, and = 

(271(7^)“^/^ exp(—C^/2cr^) = Pz[C] be the Gaussian pdf of 
noise, the conditional mean becomes 

<^[y] =2/log I cr^)I (38) 

and it reduces to (j)[y\ = tanh[y/a^] for x € {—1, +1} (BPSK 
constellation). 
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